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ABSTRACT 

A simple computational procedure has been developed by using Monte Carlo for allocating redundancy among 
subsystems so as to achieve maximum reliability of a complex systems subjected to multiple constraints which may be 
linear, non linear separable or non separable. Two examples of linear, non linear separable and non separable constraints 
with having twenty problems are solved. 

KEYWORDS: Active Redundancy, System Reliability, Constraints, Problem Formulation, Monte Carlo Algorithm 
INTRODUCTION 

Reliability of an overall system can be increased by introducing redundancies in subsystems. In order to ensure 
that factors such as cost, weight, and volume remain within resources available, system reliability is optimized with respect 
to these constraints. 

The redundancy optimization problem is a classical problem that has attracted considerable attention from the 
research community. It has been solved by using optimization approaches and techniques such as dynamic 
programming[l], integer programming[3], mixed-integer linear programming [2], heuristics[4] and meta-heuristics[5] 
which gives either optimal or very good solutions. 

This paper presents a Monte Carlo Method [6] that solves the problems of constrained redundancy optimization in 
complex system which gives the local optimal system reliability. Our aim is to show how to maximize system reliability 
subjected to the given constraints. For optimizing the system reliability we prefer simulation because it gives the result 
fast & accurate. 

For solving this problem (or optimizing system reliability) many authors have used cumbersome approaches. 
However in this paper the proposed Monte Carlo Method is applied for all systems and to all type of constraints and also it 
has been tried on many problems, with satisfactory results. 

PROBLEM FORMULATION AND COMPUTATION PROCEDURE 
Notations 

• Ri(Xj), Qi(Xi) Reliability, unreliability of subsystem-i with x ; components. 

• Rs(x) System reliability 

• Xi Number of components in subsystem-i 
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• x (x b x n ) 

• x*(r u ) Optimum value of x in iteration u 

• G ir (Xi) Resource-r consumed in subsystem-i with x t components. 

• b r Maximum of resource-r 

• k Number of constraints 

• n Number of subsystems 

• f(.) Function that yields the system reliability, based on n unique, subsystems and which depends on 

the configuration of the subsystems. 

• d] (x) Sensitivity-factor 

• bi(Xj) Selection-factor 

• pi Minimal path Set-1 of the system 

• q Cost of the X; subsystem 
Assumptions 



• The system and all of its subsystems are coherent. 

• There are n subsystems in the system. Subsystem structure (other than coherence) is not restricted. 

• All component states are mutually statistically in dependent. 

• All constraints are separable and additive among components. Each constraint is an increasing function of x ; for 
each subsystem. 

• Redundant components can not cross subsystem boundaries. 
Problem Formulation 

• The problem of constrained redundancy optimization in complex systems can be reduced to the following integer 
programming problem: 

• Maximize: Rs(x) = f (RiCxj), ,R„(x n )) 

• Subject to multiple linear or nonlinear constraints : 

2f =1 _ G-:,{x) < b T ; r = 1,2,3 s where the x- are positive integers, for all i 

MONTE CARLO ALGORITHM 
Steps of Algorithm 

Let X; =1 for all i 

Calculate Rs(x) 
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Repeat 100 times 

*i = random number between 1 to 10 for all i. 

Both linear and nonlinear (^■i=i.'^i cC-^O — &r) are violated. If no constraint are violated 
Calculate Rs(x*) 
IfRs(x*)>Rs(x) 

Then Rs(x) = Rs(x*) at x = x* 

Calculate the output G(*) and RsOO 
The simulation is limited to 3000 iteration only 
FIRST NUMERICAL EXAMPLE [7] 



r— Ri 



P, — 1 



Figure 1: A Bridge System 

The expression of system reliability (R s ) is 

R s (x) = R 1 (x 1 )R 2 (x 2 )Q 3 (x 3 )Q 5 (x 5 ) + Q,(x,) R 3 (x 3 ) R 4 (x 4 ) Q 5 (x 5 ) + [R l (x l )R 3 (x 3 ) + R 3 (x 3 )R 5 (x 5 ) + R 5 (x 5 )R 1 (x 1 ) 
2R 1 (x 1 )R 3 (x 3 )R 5 (x 5 )] *[R 2 (x 2 ) + R 4 (x 4 ) - R 2 (x 2 )R 4 (x 4 )] 



Problem 



Table 1: Shows Subsystem's Data 



i 


1 


2 


3 


4 


5 


R 


0.70 


0.85 


0.75 


0.80 


0.90 


Qi 


0.30 


0.15 


0.25 


0.20 


0.10 


Ci 


2 


3 


2 


3 


1 



Solution 



Maximize R s (x) 

Subjected to linear constraints H f=i CiX, < 20 Where Xi_ i= 0, ,5 are positive integers 

Table 2: Shows the Optimal Solution after Particular Iteration 



Ri 


0.70 


0.85 


0.75 


0.80 


0.90 


ZCiXi 


R s (x) 


Iteration (u) 


Q, 


0.30 


0.15 


0.25 


0.20 


0.10 








c, 


2 


3 


2 


3 


1 


20 


0.9949 


2400 


Xi 


1 


2 


4 


1 


1 


CjXj 


2 


6 


8 


3 


1 
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The above table shows the result of Monte Carlo Simulation. At 2400 iteration we get best combination of 

redundancy (x,, x 5 ), which satisfy the given constraints, It means the optimal solution is x*= {1,2,4,1,1} and 

R s (x*)= 0.9949. 

TWENTY DIFFERENT SETS OF PROBLEM 

Twenty different data sets are considered and are solved by the proposed algorithm to obtain the best combination 
of redundancy (Xj) and system reliability (R s ). 

Problems 

Table 3: Shows the Cost of the Subsystems 



Cost (d) 


Cl 


c 2 


c 3 


c 4 


c s 




2 


3 


2 


3 


1 



Table 4: Shows the Subsystem's Reliability and Consumed Resources 



Problem No. 


Ri(xi) 


R 2 (x 2 ) 


R 3 (x 3 ) 


R 4 (x 4 ) 


R 5 (x 5 ) 




1 


0.70 


0.85 


0.75 


0.80 


0.90 


20 


2 


0.60 


0.70 


0.80 


0.90 


0.60 


30 


3 


0.70 


0.85 


0.90 


0.80 


0.75 


30 


4 


0.70 


0.85 


0.90 


0.80 


0.75 


35 


5 


0.80 


0.80 


0.90 


0.90 


0.90 


25 


6 


0.80 


0.80 


0.90 


0.90 


0.90 


30 


7 


0.80 


0.80 


0.90 


0.90 


0.90 


35 


8 


0.75 


0.65 


0.85 


0.95 


0.55 


20 


9 


0.75 


0.65 


0.85 


0.95 


0.55 


25 


10 


0.75 


0.65 


0.85 


0.95 


0.55 


30 


11 


0.80 


0.65 


0.72 


0.78 


0.60 


25 


12 


0.80 


0.65 


0.72 


0.78 


0.60 


30 


13 


0.80 


0.65 


0.72 


0.78 


0.60 


34 


14 


0.90 


0.90 


0.75 


0.75 


0.75 


25 


15 


0.90 


0.90 


0.75 


0.80 


0.75 


25 


16 


0.90 


0.90 


0.75 


0.80 


0.75 


30 


17 


0.72 


0.62 


0.59 


0.93 


0.67 


25 


18 


0.72 


0.62 


0.59 


0.93 


0.67 


30 


19 


0.72 


0.62 


0.59 


0.93 


0.67 


35 


20 


0.95 


0.67 


0.59 


0.93 


0.67 


30 



Solutions 

Table 5: Shows Optimal Redundancy and System Reliability for the Given Problems 



Problem 


Redundancy 


Consumed 


Reliability 


Iteration 


No. 


{xj 


Resources (XcjxO 


(Rs(Xi)) 


(u) 


1 


1,2,4,1,1 


20 


0.9949 


2600 


2 


1,2,5,3,3 


30 


0.9998 


1900 


3 


4,2,1,4,2 


30 


0.9991 


2100 


4 


5,3,3,3,1 


35 


0.99996 


1700 


5 


3,1,1,4,2 


25 


0.9990 


2700 


6 


2,3,3,3,2 


30 


0.9999 


1600 


7 


4,1,4,4,4 


35 


0.99996 


800 


8 


3,1,2,2,1 


20 


0.9953 


2600 


9 


4,2,1,2,3 


25 


0.9974 


2300 


10 


2,3,5,2,1 


30 


0.9998 


1500 
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Table 5: Contd.. 



11 


3,1,2,3,3 


25 


0.9939 


2300 


12 


2,2,3,3,5 


30 


0.9978 


1100 


13 


1,2,5,5,1 


34 


0.9995 


900 


14 


3,3,2,1,1 


23 


0.9997 


2800 


15 


3,3,1,1,5 


25 


0.9995 


2100 


16 


3,3,1,4,1 


30 


0.9997 


1700 


17 


5,1,1,2,4 


25 


0.9956 


1900 


18 


3,2,4,2,4 


30 


0.9986 


1200 


19 


4,1,5,3,5 


35 


0.9998 


900 


20 


2,1,5,3,4 


30 


0.9998 


1400 



SECOND NUMERICAL EXAMPLE 




Figure 2: A Complex System (5 Nodes, 7 Edges) 



The system reliability (R s ) equation is 

R S (X) = ( Qs* s )[ (l- Ci" 1 )d- Qi" )(i- Qs" )( Q?" )( )+( Qi*> )( Q 2 " )(i- Q*~* )(i- Qb" ) 

Problem 

Maximize the System Reliability R s (x) 

Subject to the given constraints ^i=l c i x i < 35 where x i? i = 1, 7 are positive integers. 

Solution 



Table 6: Shows the Optimal Solution after Particular Iteration 



Ri 


0.60 


0.55 


0.75 


0.90 


0.90 


0.75 


0.60 




R s (x) 


Iteration(u) 


Q, 


0.40 


0.45 


0.25 


0.10 


0.10 


0.25 


0.40 








c, 


2 


3 


3 


2 


2 


1 


1 


35 


0.9899 


1700 




1 


2 


3 


3 


2 


3 


5 


CjXj 


2 


6 


9 


6 


4 


3 


5 









The above table shows the result of Monte Carlo Simulation; at 1700 iteration we get best combination of 
redundancy (X[ x 7 ). It means the optimal solution is x*= {1,2,3,3,2,3,5} and Rs(x*) = 0.9899. 

APPLICATIONS 

Monte Carlo is also applicable for both linear and nonlinear Constraints. Same bridge systems as shown in 
figure 1 are taken and new twenty different problems are framed by changing the each subsystem reliability, we can obtain 
the optimal system reliability under the subjected non linear constraints. 
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Maximize the system reliability R s (x) = R 1 (x 1 )R 2 (x2)Q3(x 3 )Q5(x 5 ) + QjCxj) R 3 (x 3 ) R 4 (x 4 ) Q 5 (x 5 ) + [R 1 (x 1 )R 3 (x 3 ) + 

R 3 (X 3 )R 5 (X 5 )+R 5 (X 5 )R 1 (X 1 )-2R I (X 1 )R 3 (X 3 )R 5 (X 5 )]*[R2(X2)+R4(X4)-R2(X 2 )R4(X4)] 

Subjected to the Non Linear (Separable) Constraints [8] 

gi(x) = x, 2 +2x 2 2 +3 x 3 2 +4x 4 2 +2x 5 2 < Ci 

g 2 (x) = 7(x 1 +exp(x 1 /4))+7(x 2 +exp(x 2 /4))+5(x 3 +exp(x 3 /4)) + 9(x 4 +exp(x4/4)+4(x5+exp(x 5 /4)) < C 2 
g 3 (x) = 7xiexp(x!/4) + 8x 2 exp(x 2 /4) +8x 3 exp(x 3 /4) + 6x 4 exp(x4/4) +9x 5 exp(x 5 /4) < C 3 
Twenty Different Problems 



Table 7: Shows 20 Problems Containing Subsystem's Reliability 



Problem No. 


Ri(xi) 


R 2 (x 2 ) 


R 3 (x 3 ) 


R 4 (x 4 ) 


Rs(xs) 


C, 


c 2 


c, 


1 


0.80 


0.85 


0.90 


0.65 


0.75 


110 


175 


200 


2 


0.70 


0.90 


0.55 


0.65 


0.85 


110 


175 


200 


3 


0.90 


0.95 


0.55 


0.85 


0.55 


110 


175 


220 


4 


0.80 


0.70 


0.60 


0.80 


0.70 


115 


140 


220 


5 


0.60 


0.78 


0.65 


0.82 


0.70 


115 


140 


220 


6 


0.60 


0.78 


0.65 


0.82 


0.70 


120 


130 


180 


7 


0.90 


0.88 


0.82 


0.83 


0.50 


120 


130 


180 


8 


0.90 


0.88 


0.82 


0.83 


0.50 


150 


140 


210 


9 


0.86 


0.79 


0.78 


0.88 


0.78 


150 


140 


210 


10 


0.65 


0.75 


0.55 


0.45 


0.85 


150 


140 


210 


11 


0.65 


0.75 


0.55 


0.45 


0.85 


125 


155 


215 


12 


0.70 


0.80 


0.90 


0.90 


0.90 


125 


155 


215 


13 


0.78 


0.88 


0.91 


0.81 


0.76 


125 


155 


215 


14 


0.78 


0.88 


0.91 


0.81 


0.76 


150 


150 


250 


15 


0.90 


0.80 


0.75 


0.65 


0.83 


150 


150 


250 


16 


0.9 


0.6 


0.76 


0.79 


0.81 


110 


175 


200 


17 


0.9 


0.6 


0.76 


0.79 


0.81 


113 


177 


210 


18 


0.77 


0.68 


0.59 


0.91 


0.86 


54 


138 


90 


19 


0.77 


0.68 


0.59 


0.91 


0.86 


61 


142 


187 


20 


0.58 


0.61 


0.93 


0.74 


0.9 


112 


147 


187 



Solutions of Twenty Different Problems 



Table 8: Shows Optimal Redundancy, System Reliability and Iteration 



Problem 
No. 


Redundancy 


Constrains 


Reliability 


Iteration 


Xl,X2,X3,X 4 ,X 5 


Ci 


c 2 


c 3 


Rs(Xi) 


u 


1 


3,3,3,2,1 


110 


175 


200 


0.9993 


2800 


2 


2,1,4,3,1 


110 


175 


200 


0.9910 


1800 


3 


4,2,2,1,3 


110 


175 


220 


0.9996 


900 


4 


3,2,2,2,2 


115 


140 


220 


0.9939 


2500 


5 


3,2,3,2,1 


115 


140 


220 


0.9933 


1900 


6 


2,2,4,2,1 


120 


130 


180 


0.9944 


2600 


7 


1,1,2,3,3 


120 


130 


180 


0.9957 


1900 


8 


1,3,3,2,2 


150 


140 


210 


0.9986 


1100 


9 


2,1,3,3,1 


150 


140 


210 


0.9989 


1700 


10 


2,4,3,1,1 


150 


140 


210 


0.9775 


1600 


11 


2,4,3,3,1 


125 


155 


215 


0.9854 


2100 


12 


3,2,3,2,1 


125 


155 


215 


0.9995 


1200 


13 


3,2,2,3,3 


125 


155 


215 


0.9998 


1000 


14 


2,3,3,3,2 


150 


150 


250 


0.9858 


1300 



Index Copernicus Value: 3.0 - Articles can be sent to editor@impactjournals.us 



Heuristic Algorithm on Monte Carlo for Constrained Redundancy Optimization of Complex System 



71 



Table 8: Contd. 



15 


3,4,3,1,1 


150 


150 


250 


0.9994 


1600 


16 


1,2,3,3,3 


110 


175 


200 


0.9971 


1100 


17 


3,2,3,4,1 


113 


177 


210 


0.9993 


1400 


18 


5,2,1,2,1 


54 


138 


90 


0.9931 


1700 


19 


4,3,1,2,2 


61 


142 


187 


0.9983 


1900 


20 


1,3,3,4,1 


112 


147 


187 


0.9994 


2400 



Subjected to Non Linear (Non Separable) Constraints 

gi(x) = 2xiX 2 2 + 3 x 3 x 4 2 + 4x 4 x 5 2 < Q 

g 2 (x) = 7(x 1 +exp(xi/4))+7(x 2 +exp(x2/4))+5(x 3 +exp(x3/4))+9(x4+exp(x 4 /4) < C 2 
g 3 (x) = 7x!exp(x 2 /4) + 8x 2 exp(x 1 /4) +8x 3 exp(x4/4) + 6x 4 exp(x 3 /4) + 9x 4 exp(x 5 /4) < C 3 
Twenty Different Problems 

Table 9: Contains the Subsystem's Reliability and Consumed Resources 



Problem 
No. 


Reliabilities R (xj) 


Constraints 


Ri 


R 2 


R 3 


R 4 


Rs 


Cx 


c 2 


c 3 


1 


0.80 


0.85 


0.90 


0.65 


0.75 


110 


175 


200 


2 


0.70 


0.90 


0.55 


0.65 


0.85 


110 


175 


200 


3 


0.90 


0.95 


0.55 


0.85 


0.55 


110 


175 


220 


4 


0.80 


0.70 


0.60 


0.80 


0.70 


115 


140 


220 


5 


0.60 


0.78 


0.65 


0.82 


0.70 


115 


140 


220 


6 


0.60 


0.78 


0.65 


0.82 


0.70 


120 


130 


180 


7 


0.90 


0.88 


0.82 


0.83 


0.50 


120 


130 


180 


8 


0.90 


0.88 


0.82 


0.83 


0.50 


150 


140 


210 


9 


0.86 


0.79 


0.78 


0.88 


0.78 


150 


140 


210 


10 


0.65 


0.75 


0.55 


0.45 


0.85 


150 


140 


210 


11 


0.65 


0.75 


0.55 


0.45 


0.85 


125 


155 


215 


12 


0.70 


0.80 


0.90 


0.90 


0.90 


125 


155 


215 


13 


0.78 


0.88 


0.91 


0.81 


0.76 


130 


155 


215 


14 


0.78 


0.88 


0.91 


0.81 


0.76 


150 


150 


250 


15 


0.90 


0.80 


0.75 


0.65 


0.83 


150 


150 


250 


16 


0.90 


0.60 


0.76 


0.79 


0.81 


110 


175 


200 


17 


0.90 


0.60 


0.76 


0.79 


0.81 


113 


177 


210 


18 


0.77 


0.68 


0.59 


0.91 


0.86 


70 


145 


150 


19 


0.77 


0.68 


0.59 


0.91 


0.86 


70 


138 


140 


20 


0.58 


0.61 


0.93 


0.74 


0.90 


170 


160 


210 



Solution of Twenty Different Problems 



Table 10: Shows Optimal Redundancy and System Reliability and Constrains 



Problem 
No. 


Redundancy 


Constraints 


Reliability 


Iteration 


{x,} 


C! 


c 2 


c 3 


Rs(xO 


u 


1 


1,3,4,1,3 


66 


126 


122 


0.9977 


1600 


2 


5,2,4,1,2 


68 


154 


186 


0.9964 


2600 


3 


3,2,2,1,2 


46 


115 


114 


0.9993 


600 


4 


2,3,3,2,2 


104 


134 


164 


0.9961 


2400 


5 


3,3,2,2,1 


86 


132 


165 


0.9909 


2700 


6 


1,2,4,2,2 


88 


123 


147 


0.9913 


1800 


7 


2,3,1,2,2 


80 


120 


127 


0.9980 


2800 


8 


3,3,2,2,2 


110 


137 


171 


0.9999 


1400 
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Table 10: Contd. 



9 


1,3,3,3,1 


111 


133 


169 


0.9984 


1700 


10 


3,3,4,1,1 


70 


135 


164 


0.9863 


2500 


11 


3,4,3,1,1 


109 


138 


180 


0.9907 


2800 


12 


4,3,1,2,1 


116 


142 


183 


0.9991 


1400 


13 


2,3,3,3,1 


129 


142 


193 


0.9999 


1800 


14 


1,3,3,3,2 


147 


138 


179 


0.9997 


1900 


15 


2,3,3,3,1 


129 


142 


193 


0.9994 


2400 


16 


3,3,3,2,1 


98 


139 


183 


0.9970 


1200 


17 


1,3,3,3,1 


111 


133 


169 


0.9977 


2600 


18 


5,1,2,2,2 


66 


141 


149 


0.9963 


2100 


19 


1,3,5,1,3 


69 


135 


137 


0.9943 


2800 


20 


4,3,2,2,3 


168 


168 


209 


0.9992 


1800 



RESULTS AND CONCLUSIONS 

In this paper, the Proposed Monte Carlo method is applied to different problems of bridge and complex 
(5 nodes, 7edges) system to optimize the system reliability under the given constraints. 

This leads to the conclusion that the system reliability of Monte Carlo is much more optimized and we can get 
better result. The proposed method is also not affected by any kind or number of constraints which is shown by applying 
different problems. Its simplicity and suitability for computer solution makes this method highly useful to reliability 
engineers. It has not been rigorously proved that this method is exact, but no faulty solution has yet been found. Even if the 
method is not exact, it appears to be very close and we get better optimal solution for any complex system. 
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